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Abstract 

We derive an exact, analytic expression for the fourth virial coefficient of a 
system of polydisperse spheres under the constraint that the smallest sphere 
has a radius smaller than a given function of the radii of the three remaining 
particles. 

I. INTRODUCTION 

It is surprisingly difficult to calculate analytic expressions for virial coefficients. Even the 
exact calculation of the second virial coefficient is in general extremely difficult. If particles 
have no other interaction than hard core repulsion, things are slightly easier. In the isotropic 
phase, the second virial coefficient B2 of two arbitrary, convex particles, A and 5, can be 
determined by [Q 

B2{A, B) = ^ {y{A)G{B) + S{A)M{B) + M{A)S{B) + G{A)V{B)) (1) 

OTT 

where only geometrical quantities of both particles are used: the volume V , the surface area 
S, and the integrals M and G over the mean respectively Gaussian curvature. For higher 
virial coefficients only approximations of a similar type are known Q]. 

The second and third virial coefficients of identical spheres are known analytically and 
are obtained after a straightforward calculation. The fourth virial coefficient, for which the 
analytic expression due to Boltzmann is also known 0, is much more difficult. Beyond 
the fourth virial coefficient only numerical data are available Q. More recently, data of the 
fourth and fifth virial coefficients of binary mixtures have been published by Saija et al 

In this paper we will focus on hard spheres with different radii, for which until now, 
no analytic results are known beyond the third virial coefficient . In section [I^ we show the 



results for the second and third virial coefficients. In section III we derive the known results 



for the fourth virial coefficient. The main problem, the calculation of the complete star 



diagram, is derived in section together with an inequality, which has to be satisfied, in 
order for the given expression to be valid. In section |V| we finish by discussing our result 
and some concluding remarks. 
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II. SECOND AND THIRD VIRIAL COEFFICIENT 



The general expression for the nth virial coefficient of a gas with pairwise, additive 
interaction (pjj between particles i and j is given in terms of Mayer f^j functions by 



B. 



lim ¥-'[■■■[ Vndr,---drn, (2) 

^^oo J J 



n! V 

where r.; arc the spatial coordinates of particle i, and V is the volume accessible for the 
particles. is the sum over all labeled stars with n points given by 

n 

K ^ E n k (3) 

{Sn} i>3 

and expressed in the Mayer / functions, which are related to the interaction potential by 

fij = exp(-/50ij) - 1. (4) 

For hard particles the Mayer / function reduces to 

f — ] if no overlap , 
1^ —1 if overlap ^ ' 

because the interaction potential is either zero or infinity in the case that particles do overlap, 
respectively do not overlap. 

The second virial coefficient B2{A, B) for two particles A and B, has only one contributing 
diagram, and is given by 

S2(A5) = ^ OO =^ , (6) 

where we have explicitly labeled the diagram. We use diagrams to represent the different 
integrals. Integrations are over all possible positions of the different particles with one parti- 
cle fixed. Each solid line represents an overlap between the connected particles, and , hence, 
some constraint on the position vectors. The signs, which are determined by the Mayer 
functions, are, for convenience, immediately put in front of the diagrams. The diagram in 
B2 is the volume from which the second particle is excluded, in order not to cause an overlap. 
For spheres with radii A and B this is simply the volume of a sphere with radius A-\- B, 
leading to the well known result 

B2{A,B) = %,{A + Bf. (7) 

Also the third virial coefficient B^{A,B,C) of three particles A, B and C has only one 
contributing diagram 

BM,B,C) = l ^ =\ ^. (8) 
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Only in the case of simple objects, e.g. spheres and discs, are analytic expressions for the 
third virial coefficient known. 

In order to proceed we will assume that the radii of the three spheres A, B and C 
satisfy the inequality A > B > C . We position the biggest particle at the origin and 
particle B on the positive z-axis. The integral corresponding to the diagram consists of 
two contributions. The first contribution comes from the case where B is completely inside 
particle A, in which case particle C is only required to overlap with particle 5, leading 
automatically to an overlap with particle A. The second contribution comes from the case 
where B is overlapping with A, but not enclosed. Particle C now has to overlap with both 
particles, A and -B, and hence its center has to be placed in the overlapping volume of the 
spheres placed at A with radius A + C and the sphere placed at B with radius B + C. This 
volume Z{zab, A + C,B + C) is the sum of two segments of different spheres and depends 
on the distance za_b between A and B. The formula Z{r, Ri, R2) can easily be evaluated if 
we assume Ri > R2 

r fRl r<Ri-R2 
Z{r, R,, R2) = \ ^^{Ri + R2- r)2(r2 - 3(i?i - i?2)' + 2r(i?i + R2)) (9) 
[0 r> R1 + R2 

if r < Ri — R2 the smaller sphere is completely enclosed by the larger sphere, and, hence, 
the overlap is equal to the volume of the smaller sphere. If r > Ri + R2 the spheres are not 
overlapping at all. The diagram for the third virial coefficient is therefore given by 

A+B 

= J dzAB^T^Z%Z{zAB,A + C,B + C) (10) 


which results for the third virial coefficient B^{A, B, C) in 

B^{A, B, C) = -^{A^B^ + B'^C^ + A^C'^ + 

^ABC[A + B + C]\AB + BC + CA]). (11) 

Although we have assumed that A > B > C, the resulting formula is fully symmetric in the 
radii. 



III. FOURTH VIRIAL COEFFICIENT 

In order to obtain the fourth virial coefficient B4 for particles A, B, C and D, we need 
to evaluate three different diagrams, although this could be reduced to only two modified 
star diagrams (see ref. ^j) 



B,{A,B,C,D) = - 




The fourth virial coefficient is here expressed in unlabeled diagrams. However, labeling now 
becomes important because we want to evaluate those diagrams for the case of non-identical 
particles. This leads to the following expressions 
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(d ) — ( c ) f b) — ( d) (c) — f b) ( d) — f c) (a) — ( d) f b) — ( d) 

( A ) ( B ) ( A ) J ( A ) ( D J ( B ) ( A ) ( B ) ( C J ( C ) ( A ) 





In order to proceed, we need to assume that the radii of the spheres satisfy A > B > 
C > D. The labeled diagrams in the first line of (|13]) are given by expressions of the form 



A+C+2D 

(D ) (C ■ 





j dzAc^7rz\cZ{zAc,A + B,C + B)Z{zAc,A + D,C + D). (14) 



The limits of integration are a consequence of the differences in size and in this case deter- 
mined by the fact that B > D. The diagrams in the second line of (|1^) are similar, and 
only differ in the integration limits because of the extra overlap. They are of the form 

A+C 

J dzAc'^7rz\cZ{zAc,A + B,C + B)Z{zAc,A + D,C + D). (15) 



The other diagrams are obtained by permutation symmetry. The resulting integrals of ([M]) 
and ([15|) can easily be evaluated, but lead to lengthy expressions. The summations over 
the different labelings do not lead to symmetric expressions in terms of the radii. If we 
however take the combination of the diagrams according to the definition of the fourth virial 
coefficient (|12|), we obtain 

6 1/X = {A^B^C^ + A^B^D^ + A^C^D^ + B^C^D^ + 

O — O 9 

3ABCD[AB + AC + BC + AD + BD + CD] [ABC + ABD + ACD + BCD]), (16) 
which is symmetric in the different radii of the spheres. 

IV. THE COMPLETE STAR DIAGRAM 

The problem of calculating the fourth virial coefficient lies in the remaining diagram, the 
complete star diagram in which all pairs of the four particles overlap. An analytic expression 
for this diagram in general is not known. What is known is the special case of four identical 
radii, which is due to Boltzmann 

2^67r^i?^ 

( 34197r - 438^2 - 8262 arccos(l /v^) ) , (17) 

945 V V / V ;y , y J 
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where we used R to denote the radius of the particles. 

Some hmiting cases for different radii are easily obtained. In the limit that the radius of 
the smallest sphere goes to zero the three remaining particles have to overlap with a point 



lim 




(18) 



In the limit where the radius of the biggest sphere goes to infinity, the other three particles 
give rise to the diagram related to 



lim 




A' 



(19) 



In the limit where both A and B go to infinity 



lim 

A.B^oo 




lim 

A,B-*oo 



A^B^{C + Df 



(20) 



In order to explore the star diagram for non-limiting cases we must look in more detail to 
the diagrams and introduce the /-bonds as defined by Ree and Hoover ^ which we denote 
by a dashed line. Such a connection refers to non-overlapping particles, solid connections to 
overlapping particles and no connection allows for both, overlap and no overlap. 

We will now consider diagrams in which the particles A, B and C all overlap with the 
smallest particle D. The diagram without other constraints can be written as the following 
summation of diagrams: 





+ 




+ 




+ 




+ 




(D) (C) (D) (C) 

+ 






(21) 



where we have denoted at the right hand side explicitly whether pairs of particles overlap 
or not. The diagram at the left hand side can, however, be evaluated immediately 




(22) 



The same is true for combinations of the form: 






( D ) ( C ) ( D } ( C } (p ) (CJ (p J (C) [pj (C) 

+ JXJ. + XXX + IXT = (£M2J X 





(23) 
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As well, some combinations have already been evaluated before (|T5D, e.g. 








+ JNX + Xak - ^Xa • (25) 



For all diagrams on the right hand side we have an analytic expression with only one excep- 
tion - the last diagram. We demonstrate below that this last diagram is identical to zero if 
the radii of the spheres satisfy a simple inequality. 

In order for this diagram to be non-zero the spheres with radius A, B and C should all 
overlap with the smallest sphere with radius but have no overlaps with each other. If the 
spheres A, B and C are mutually non-overlapping as shown in figure |I], there is a smallest 
sphere D that can be constructed that will touch with all three spheres. The centers of mass 
must all lie in the same plane, and therefore this problem reduces to the one of touching 
circles, known as Apollonian circles since the construction of this fourth circle was first 
solved by the Greek mathematician ApoUonius of Perga (200 BC) . If the radius D were 
smaller, at least one of the spheres A, B oi C would have no overlap with this circle. 

If this were the case for all possible configurations of non-overlapping spheres A, B and 
C there would be no contribution to the last diagram of (0), and hence we would know the 
fourth virial coefficient. The smallest radius D which can contribute is found in the case that 
all four spheres touch as depicted in figure ^ For this special case of the Apollonian problem, 
a simple relation between the radii of the circles was derived by the French philosopher and 
mathematician Descartes [ITUl, known as the Descartes circle theorem: 



I^B + C + fl) =Hl^+B5^CJ+fl5)- (26) 



Solving this equation gives the upper limit for radius D 



O < (27) 

AB + AC + BC + 2JaBC{A + B + C) 



which will lead to a zero- valued diagram, and hence enables us to give an analytic expression 
for the complete star diagram. Therefore we obtain an analytic expression for the fourth 
virial coefficient of polydisperse spheres under the constraint (pTl) 

B^iA, B, C, D) = (i|f^) X (A^E^C^ + ^353^3 + ^3^3^3 + 53^3^3+ 

2,ABCD[AB + AC + BC + AD + BD + CD] [ABC + ABD + ACD + BCD])- 
^D'^A^B^C^ + ^D^ABC{AB + AC + BC)- . . 

^D\A^B^ -2A^BC -2AB^C + A^C^ -2ABC^ + B^C^)- ^ ' 

^D\A + B){A + C){B + C)~ ^D\A^ + ^AB + B^ + 3AC + 35C + C^)- 
^D^{A + B + C)-^D^ 
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The first part is symmetric in the four different radii and proportional to ([I6D, the last part, 



however, is only symmetric in the radii A, B and C. Note that the limiting cases (18), (19) 
and ( |20| ) are in agreement with this result. 

In the case that A = B = C = 1, and, since according to (|3) D < (273-3)73 ^ 0.1547, 
this formula reduces to 

54(1, 1, 1, D) = (g) (70 + 630D + 2520^2 + U70D^ + 945D^+ 

189D^ - 1512D6 - 1296D^ - 405D^ - 450^) . (29) 

This result can be compared with the numerical data for the binary mixture |^. Only the 
first three numbers in the first column satisfy the inequality (^) 

Dni2in = i) = S5— 000 = 0.0553630659 [0.05539(2)] 

All2(7^ = |j) = Soo—o = 0.0819065785 [0.08189(2)] (30) 

Diiuin = A) = f^SS: = 0.1161886732 [0.11617(3)] 

The numerical values obtained by Saija et al. [0 are given in square brackets. 



V. DISCUSSION 

To our knowledge this is the first time that an exact and analytic solution of the fourth 
virial coefficient of a hard sphere mixture is found. One should note, however, that in order 
for the expression to be valid, the radius of the smallest sphere is at most 0.1547 times that 
of the largest sphere. Our result agrees perfectly with the relevant, numerical results of 
Saija a/ 1]^. However, in the description of a polydisperse mixture of spheres, the most 
interesting case is that of a few large and many small particles. Unfortunately, this also 
means that the inequality (|27| ) is, in general, not satisfied. The inequality requires that 
the smallest radius is strictly smaller than the other three, which, for instance in a binary 
mixture, allows only one of the three mixed virial coefficients to be evaluated. 

The method described in this article can almost certainly be used to evaluate the fourth 
virial coefficient of other particles, e.g. discs, and can possibly be used to simplify calcula- 
tions of higher virial coefficients of asymmetric mixtures. 
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FIGURE CAPTIONS 



1. The Apollonian problem of constructing a circle D which is tangent to three arbitrary 
spheres A, B and C. 

2. A special case of the Apollonian problem in which all spheres are tangent to each other, 
and for which case a simple relation (^Bj) between the radii was derived by Descartes. 
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FIGURES 




FIG. 1. Blaak 




FIG. 2. Blaak 
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